Spin dynamics of the quasi two dimensional spin-| quantum magnet Cs 2 CuCl 4 
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We study dynamical properties of the anisotropic triangular quantum antiferromagnet CS2CUCI4. Inelastic 
neutron scattering measurements have established that the dynamical spin correlations cannot be understood 
within a linear spin wave analysis. We go beyond linear spin wave theory by taking interactions between 
magnons into account in a 1/S expansion. We determine the dynamical structure factor and carry out extensive 
comparisons with experimental data. We find that compared to linear spin wave theory a significant fraction 
of the scattering intensity is shifted to higher energies and strong scattering continua are present. However, the 
1/S expansion fails to account for the experimentally observed large quantum renormalization of the exchange 
energies. 
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I. INTRODUCTION 

The quasi two dimensional spin- 1/2 quantum magnet 
CS2CUCI4 has attracted much theoretical and experimental in- 
terest in recent years as a possible realization of a two dimen- 
sional quantum spin iiq U idiIi2iii5i!iL2i2ii2iiiii2, This anisotropic 
triangular Heisenberg antiferromagnet is believed to be a 
promising candidate due to its small spin, quasi two dimen- 
sionality and geometrically frustrated spin interactions. Al- 
though CS2CUCI4 exhibits conventional incommensurate long 
range magnetic order at low temperatures, neutron scattering 
measurements have revealed unusual features in the spin exci- 
tation spectrum. In particular, the dynamical correlations are 
found to be dominated by an extended scattering continuum 
over a relatively large window of energies. Several workers 
have interpreted this observation as a signature of deconfined, 
fractionalized spin- 1/2 (spinon) excitations, characteristic of a 
spin liquid phase. In this line of approach, the observed broad 
scattering continuum is interpreted in terms of a two-spinon 
scattering continuumii2ii£ 

However, a strong scattering continuum does not entail an 
underlying spin liquid phase. In fact, a conventional mag- 
netically ordered phase with strong magnon interactions can 
exhibit a broad continuum due to multi magnon scattering pro- 
cesses. A previous examination of the inelastic neutron scat- 
tering data on CS2CUCI4 was performed in the framework of 
linear spin wave (LSW) theory. 4 The latter predicts sharp sin- 
gle particle excitations and weak two magnon scattering con- 
tinua, features which were argued to be in poor agreement 
with the data. Given that the magnetic properties derive from 
small S — 1/2 Cu spins, one would a priori expect magnon 
interactions to play an important role. In order to assess the 
applicability of a spin wave based scenario to CS2CUCI4 it is 
therefore necessary to go beyond linear spin wave theory. 

On a qualitative level the predictions of nonlinear spin wave 
theory are readily anticipated. By Goldstone's theorem the 
breaking of a continuous symmetry in a magnetically ordered 
state enforces the presence of single particle excitations at low 
energies. As a result of the aforementioned interactions, these 



magnons acquire a finite life time, which in turn leads to a fi- 
nite line width in the dynamical structure factor. Furthermore, 
compared to linear spin wave theory, spectral weight is trans- 
ferred to higher energies via multi magnon scattering pro- 
cesses. In the case of CS2CUCI4 one may expect the presence 
of a strong scattering continuum in the ordered phase because 
( 1 ) the low spin and the frustrated nature of the exchange in- 
teractions lead to a small ordered moment and strong quantum 
fluctuations around the ordered state; (2) the magnon interac- 
tions in non-collinear spin structures like the ones found in 
CS2CUCI4 induce a coupling between transverse and longitu- 
dinal spin fluctuations. This interaction provides an additional 
mechanism for damping the spin waves and can enhance the 
strength of the scattering continuum. 

There is evidence of low-energy spin wave modes in the in- 
elastic neutron scattering data. Sharp peaks are also observed 
at high energies near special wave vectors where a putative 
spin wave dispersion is at a saddle-point. It is important to 
note that this spin wave dispersion is dramatically "renormal- 
ized" compared to the prediction of linear spin wave theory. 3 ' 4 

A priori it appears that nonlinear spin wave theory could 
have the necessary ingredients to account for the spin corre- 
lations observed in CS2CUCI4. The issue then is whether it 
is possible to achieve a quantitative description of the exper- 
iments in low orders of perturbation theory in the spin wave 
interactions. 

In the present work we go beyond linear spin wave theory 
and include, within the framework of a 1/5 expansion, the 
quantum fluctuations around the classical ground state. We 
then apply the results to the case S = 1/2, in which the 
formal expansion parameter becomes of order 1 and is there- 
fore not small. We are motivated by the observation that spin 
wave theory gives a good description of physical properties of 
the square-lattice spin-i Heisenberg Hamiltonianpi^iii^ In- 
deed, higher order (in a 1/S expansion) corrections to linear 
spin wave theory were shown to be small in this case. Fur- 
thermore, taking these corrections into account in the calcula- 
tion of static and dynamical properties leads to an improved 
agreement with the results of more sophisticated numerical 
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techniques. l&iLAs Although a corresponding analysis is not 
available for the frustrated triangular antiferromagnet, pertur- 
bative expansions in 1/S have shown the renormalization due 
to quantum effects is relatively small>I2i22i2ii2£ 

This paper is organized as follows. The spin Hamiltonian 
for CS2G1CI4 is introduced in Sec.|ll] In Sec. Mil we determine 
the magnon Green's function in the framework of a large-S ex- 
pansion. In Sec. lIVI we relate the experimentally measured dy- 
namical correlation functions to the magnon Green's function. 
The results of our analysis and comparisons to the experimen- 
tal data on Cs 2 CuCi4 are presented in Sec.|V] We conclude 
with a summary of our results in Sec. lVIl 

II. SPIN MODEL 

The full spin Hamiltonian of CS2CUCI4 has been deter- 
mined previously from measurements in high magnetic fields 



Here the vectors Si and 5? connecting neighboring sites are 
shown in Fig.^ The vector D = (D, 0, 0) is associated with 
the oriented bond between the two coupled spins connected 
by <5i or 62 and n is a layer index. The factor (—1)™ indicates 
that the interaction alternates between even and odd layers, 
which as a result can be considered to be inverted versions of 
one another. A weak interlayer interaction J" is also present 
between neighboring layers. However, as J" is quite small we 
neglect it in the following. 
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FIG. 1: (Color Online) The magnetic sites and exchange couplings 
within a single layer of CS2C11CI4. Layers are stacked along the crys- 
tallographic a-direction with an interlayer spacing a/2 and a relative 
displacement in the c-direction. 



(see Ref. for details). For our purposes it suffices to note 
that the magnetic Cu 2+ ions form a triangular lattice with 
anisotropic exchange interactions. As shown in Fig. ^ the 
main exchange interaction J = 0.374(5) meV is along the 
crystallographic b axis ("chain direction"). A weaker spin ex- 
change J' = 0.128(5) meV occurs along the zig-zag bonds. 
Finally, a Dzyaloshinskii-Moriya (DM) interaction22i2i D = 
0.020(2) meV is present along the zig-zag bonds. 



Denoting the spin-i operators at the sites R by Sr, the 
quasi two dimensional Hamiltonian takes the form 



Sr+* 2 ) - (-l) n D • S R x (S R+il + S R+ i 2 ) . (1) 
I 




FIG. 2: (Color Online) The reciprocal space diagram of CS2C11CI4 
projected along the (0, k, I) plane. The F points refer to the center 
of the Brillouin zone and Q is the ordering wave vector. The path of 
the cut shown in Fig.|3|is depicted as a dashed line. 

Following the conventions of Coldea et al. in Ref. 0, we 
will discuss the dynamic response in terms of the two dimen- 
sional Brillouin zone of the triangular lattice even though the 
full crystal symmetry is orthorhombic. In our notation wave 
vectors are expressed in terms of the reciprocal lattice vectors 
as k = (h, k, I), which is a shorthand for 2ir(h/a, k/b, l/c). 

The Fourier transforms of the exchange and DM interac- 
tions are 

Jq = Jcos(2-7rfc) + 2 J' cos(-7rfc) cos(7rZ), (2) 
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and 



Dq = — 2iD sin(7rfc) cos(trl). 



It is convenient for what follows to define a quantity 
J Q = J Q- iD Q- 



(3) 



(4) 



Experimentally, spiral magnetic long range order is observed 
in CS2CUCI4 at temperatures below Tjy = 0.62(1)^. The 
ordered structure is found to lie in the be plane by virtue 
of the small easy-plane anisotropy generated by the DM in- 
teractions. The spin structure is an incommensurate cycloid 
with an ordering wave vector Q = (0.0, 0.5 + e, 0) where 
e = 0.030(2). 



two groups, differing according to the direction of the DM 
vector. For the case where the layer index n is odd (even), the 
DM vector is taken to point into (out of) the be plane. 

In what follows we present the results for the even layers 
only. However, it is easy to see that the spin structure factor is 
in fact independent of the layer index and the overall result is 
a simple summation over all layers. 

It is convenient to define a local reference frame (x,y, z) 
such that the classical spin direction is aligned along the z 
axis at every site 



'10 
cos(Q • R) - sin(Q • R) 
>0 sin(Q-R) cos(Q-R) 



The Holstein-Primakoff transformation reads 25 




(5) 



HI. LARGE S EXPANSION 

We now turn to a summary of our calculations. The pro- 
cedure we follow is standard. We first express the fluc- 
tuations around the "classical" ground state in terms of 
boson operators using the Holstein-Primakoff transforma- 
tion i^SiSiiSii 2 ^ 2 ^ 2 ^ 2 ^^ The term quadratic in the boson opera- 
tors constitutes the basis for linear spin wave theory, whereas 
higher order terms represent spin wave interactions. The in- 
teraction vertices of n bosons carry a factor 5 2 ~™/ 2 , where 
S is the "length" of the spin. In the second step we deter- 
mine the renormalized magnon Green's function by calculat- 
ing the self-energy to leading order in 1/5. Finally, the ex- 
perimentally observable dynamical correlation functions are 
expressed in terms of the Green's function of the Holstein- 
Primakoff bosons. 

The classical ground state is determined by treating the 
spins as classical vectors and then minimizing the energy. In 
this way one obtains a cycloidal structure with a characteristic 
wave vector Q that is fixed by the condition that it minimizes 
the exchange energy per spin, i.e. Jq = min q Jq . We find 
Q = (0.0,0.5 + e ,0) with e = 0.054. This value differs 
significantly from the measured incommensuration but quan- 
tum fluctuations lead to a reduction in eo and taking them into 
account yields good agreement with experiments^^ 

As we have already indicated in Eq. Q, to a good approxi- 
mation the layers are decoupled. Hence we consider from now 
on a set of independent 2-D layers, which are subdivided into 
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(6) 



where the boson creation and annihilation operators satisfy the 



canonical commutation relation 



'R' VR> 



Sr.r'- Here 9 



is an arbitrary angle which we set equal to jt/2 in order to 
make contact with the notation used in Ref. yl[ Introducing 
the Fourier transform 



^ R 



ikR 



(7) 



on a lattice of N sites, the Hamiltonian of Eq. Q takes the 
form 



(8) 



where Ti. n is proportional to S 2 ~ n l 2 and consists of normal 
ordered products of n boson operators. There is no Hi term, 
because Eq. (|8jl is an expansion around a minimum of the clas- 
sical energy. Linear spin wave theory takes into account only 
the terms Ho and H-2- The higher order terms represent inter- 
actions between magnons. The leading terms in the expansion 
are 
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Ho 

h 3 



NS 2 JT, 

NSJ^ + S^Ak (4 
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I 



Here the sum over k is performed in the first Brillouin zone 
and the subscripts 1 ... 4 denote It! . . . k 4 . The quantities 
Ao Bk and Ck are expressed as 



.4 k 
Bk 



J (2 J k 

J Q+k 



J Q+k 



J Q+k 
J Q-k- 



^Q-k) 



•^Q-k) 



(13) 



The coefficients and By are even functions of k, whereas 
Ck is an odd function of k. In the absence of easy-plane 
anisotropics, i.e when D vanishes and inversion symmetry is 
present, we recover the results of Ref.^J [Note that our def- 
initions in Eqs. Jl 31 differ from those of Ref. 21 by a factor 
of four.] We emphasize that the cubic interaction is generated 
as a result of the coupling between transverse and longitu- 
dinal fluctuations and hence can only exist in non-collinear 
spin structures. Furthermore, we note that the vertex factor 
Ck oc |k| 3 for small k owing to the fact that Jq is at a mini- 
mum by construction. 

The quadratic Hamiltonian H% is diagonalized by a Bogoli- 
ubov transformation 



where 



0k = w k 7 k + u k 7L k , 



A 



(14) 



The diagonal form of the quadratic Hamiltonian is 



H 2 



NSJ^ 



E 



I 7 k 7 k + § 



(15) 



(16) 



where Wk = 2SyJ A k — _B k is the linear spin wave dispersion 
relationi^Si We note that Uk is an even function of k, de- 
spite the absence of inversion symmetry in the Hamiltonian. 



In fact, the symmetry of Wk is a consequence of time-reversal 
symmetry, which implies the following relation between the 
elements of the dynamical structure factor (Eq.l27t>. 31 



°k.w 



S" 



(17) 



The importance of quantum fluctuations can be gauged by de- 
termining the average value of the local spin given by the stan- 
dard formula 



(S& = S - AS = S - ±J2 u t 

k 



(18) 



The boson Green's function at zero temperature is expressed 
as 



f>OQ 

-/ I die* 



T 



M 

-k(i) 



0k(O)0- k (O) 



(19) 



where T denotes time ordering and (...) represents a ground 
state expectation value. The inverse of the unperturbed 
Green's function is given by a 2 x 2 matrix, 



G 



(0)-l 



(-2SA k + irj)a a + 2SB k a x + luct 2 



(20) 



Here a and er denote the identity and Pauli matrices respec- 
tively and i) = + . 

The self-energy is defined by the Dyson equation, 



G 



(0)-l 



J k,w 



and can be parameterized as 

Sk,w = OkMff + X kuJ <T x + Z klAj a z 



(21) 



(22) 



The leading order (in 1/S) contributions to the self-energy can 
be divided into two parts 



,(3) 



(23) 



Here S k denotes the vacuum polarization contribution 
that arises in first order perturbation theory in H4 . It is fre- 

(3) 

quency independent and purely real. On the other hand, E k ' 
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denotes the contribution in second order perturbation theory 
of the three-magnon interaction 7i 3 . It incorporates the ef- 
fects of magnon decay. Using Eq. J20l >. the contribution 
to the self-energy is found to be of the form 



O 



(4) 



X, 



(4) 



,(4) 



2S^ 1 



l 



k' 

(A k _ k , - 5 k _ k , - A k , - A k ) A k 
25 



k' 



1 



^4k-k' — -B k _ k ' — Aw — -A k ) By 



(24) 



The contribution E^ 3 ) is most easily evaluated in the Bogoliubov basis (7) and is equal to 



O 



(3) 
k,u; 

(3) 
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1 



1 



Wk' + Wk-k' — UJ — IT) CJ k ' + ^k-k' +U) — IT) 
21 / 1 1 



-5' 



16iV 



2 {2$W (k ' j k _ k ')$(2) (k ' k _ k ')| 



^k' + ^k-k' - UJ ~ IT] U) k > + W k -k' + UJ — IT] 
1 1 



uJw + w k - k ' — lu — ir\ Uk> + w k - k ' + ui — irj 



(25) 



where 



<J> (1) (k',k - k') = (C k ' + C k _ k /) (u k < + vw) (uk-k' + «k-k') - 2C k (u k 'U k -k' + ^k'Wk-k') , 
$ (2) (k',k - k') = Ck' (uk' + «k') (wk-k' - i^-k') + Ck-k' (u k - k ' + u k -k' ) (W - w) . (26) 
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IV. DYNAMICAL CORRELATION FUNCTION 



Inelastic neutron scattering experiments probe the dynam- 
ical structure factor Sj^. The latter is defined as the Fourier 
transform of the dynamical spin-spin correlation function 



f°° elf 

S»l= e -^(S\(0)S£(t)). (27) 

J-00 27T 

Here /i, v = (a, b, c) label the various crystallographic axes 
and the Fourier-transformed spin operators are defined by 

°k - 7^ 2^r 

It is convenient to introduce time-ordered spin-spin corre- 
lation functions in the rotated coordinate system 



dte 



! <T^ k (0)^(i)>, 



(28) 



where a, (3 = (x, y, z) are the rotated coordinate axes (Eq.[5}- 
The dynamical structure factor is related to the imaginary part 



of the time ordered correlation function in the following way 

1. 



S, 



1. 



'k,u 



-Im 
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k+Q,w 



e 



k-Q,w 



nbc 

J k.w — 



= Im 

7T 



M k+Q,w W k-Q,w 



where 



©it 



pzz , z^y , ■ / pzy _ pyz \ 



(29) 
(30) 
(31) 

(32) 



To proceed further, we expand the dynamical correlation 
functions in inverse powers of S to order 0(5°). The cor- 
responding results have been derived previously by Ohyama 
and Shiba. 21 Here we merely quote their results for the sake 
of completeness. The transverse correlations are 



K% = | C 2 Tr[(a°-^)Gk,.] 



|c 2 Tr[(a° + a*) G k , u ], 



(33) 
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where the Green's function is given by Eq. (12 U and where 



^K 2 



(34) 



must be retained. The mixing of transverse and longitudinal 
fluctuations manifests itself in 



{^>[(-°+- x )G k ,.] 
+ P£Tr[^G k , w ]}. (35) 



We note that when squaring (I34> only terms to order 0(5 x ) Here the functions P^'^ are defined as 



p(i) _ _s_ 

k,u> 4JV 



^ (k',k - k') (tt k /u k _ k / + Wk»t*k-kO 

k' 

p k!i = w E $(2) ( k '> k - k ') ("k'Vk-k' + ivuk-kO 



i 



<*V + Wk-k' — uj — iri lu^i + cJ k _ k < + lu — iri 
1 1 



F, 



(0)zz _ 



-gTf E ("k'^k-k' + ^k'"k-k') 2 



1 



Finally the longitudinal correlations are decomposed in inverse powers of 5 as F k ^ = F^^ z + F^^ z , where 

1 



Wk' + w k _ k ' — uj — vq Wk' + w k _ k ' + w — 177 



Wk' + Wk-k' — u> —vq uj k / + w k _ k ' + w — ir] 

2 



= ^ { (42) Tr [(o* + o-) G k)W ] + (Pg) Tr [(a - o*) G k ,.] + 2P£p£lr [cr z G kj[ 



r 



(36) 



(37) 
(38) 



We note that the F^ zz term does not require the knowledge 
of the bosonic self-energy and is basically a free boson result. 
For this reason, it is often included in linear spin wave calcu- 
lation as a source of two magnon scattering, even though it is 
formally a higher order contribution in 1/5. In what follows, 
we abide by this (in some sense inconsistent) convention and 
consider the contribution of Eq.|^]as part of linear spin-wave 
theory. As a consequence we then retain the F^ zz contri- 
bution to the dynamical structure factor, although of higher 
order in 1/5 (i.e. 0(S~ 1 )) than the other terms we take into 
account. 

The (unpolarized) inelastic neutron scattering cross section 
is given by 

^ = l/k| 2 £(v-V^sC' 



|/ k | 2 [(i-k*)s£ 



l + K 5 



2\ ebb 



(39) 



where k M is the /i-component of the unit vector in k direction. 
The magnetic form factor / k is determined by the magnetic 
ions. For Cu 2+ , the isotropic form factor has a relatively weak 
wave vector dependence within the first Brillouin zone and 
will be neglected from now on. 32 

It is well known that the 1/5 expansion preserves many 
physical properties "order by order" in 1/5. For instance, 
it follows from Eqs. ( 1241 125> that the Goldstone modes per- 
sist beyond linear spin wave theory as one expects on physi- 
cal grounds. A careful examination also shows that to order 
0(5°) the spectral functions are positive and that the rela- 



tion illl holds. However, due to a lack of self-consistency 
the 1/5 expansion leads to an (unphysical) unequal treat- 
ment of the one-magnon and two magnon scattering contri- 
butions to dynamical correlation functions^ It is worthwhile 
to discuss this issue in more detail. The leading order con- 
tribution to the dynamical structure factor is due to coher- 
ent single magnon excitations and is of the form 6 (to — w k ). 
The two magnon contribution due to longitudinal fluctua- 
tions (Eq. I37t gives rise to a scattering continuum of the 
form / (k, k') 5 (u> — cj k < — o-> k _k') with some function 
I (k, k'). The extent of the two magnon contribution in k — oj 
space is determined by the lower and upper bounds of the 
function w k / + w k _k' for a given k. 

On general grounds, we expect the lower bound of the two 
magnon scattering continuum to be equal to or smaller than 
the "true" magnon dispersion u) k . In fact, the existence of a 
zero-momentum Goldstone mode guarantees that there exists 
a two magnon contribution at frequencies <D k + ujq = u) k . 

It is easy to see that this property does not hold order by 
order in a 1/5 expansion. Indeed, the first order contribution 
in 1/5 shifts the pole of the Green's function and leads to a 
renormalization of the magnon dispersion. The renormalized 
dispersion cD k can be determined from the Dyson equation 



- S k,^ = 0- 



(40) 



However, to order 0(5°) the threshold of the two magnon 
contribution is still determined by the bare dispersion rela- 
tion Wk- This results in an unphysical behavior, where the 
two magnon scattering continuum is separated from the single 
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magnon dispersion by a gap. In order to avoid this problem, 
we impose the following self-consistency condition: the lin- 
ear spin wave dispersion Wk used in Eqs. f25l 061071) is to be 
replaced by the renormalized dispersion u>k- 



V. DYNAMICAL PROPERTIES OF Cs 2 CuCl 4 

So far our discussion of the 1/5 expansion has been fairly 
general. In order to make contact with the experiments on 
CS2CUCI4 we now set the exchange constants to their appro- 
priate values 3 ' 4 and fix 5 = 1/2. We then evaluate the dy- 
namical structure factor at a given wave vector numerically. 
Complex integrals such as Eqs. J25II36> are evaluated by sum- 
ming the imaginary part of the integrands over a frequency 
grid of 1200 points and of 1000 x 1000 points in wave vector 
space. The real parts are then determined from the Kramers- 
Kronig relations. The aforementioned self-consistency condi- 
tion is implemented by calculating the full Green's function 
iteratively on a 100 x 100 grid in the Brillouin zone. We 
observe satisfactory numerical convergence after about three 
iterations. 

We first turn to the magnon dispersion. The linear spin 
wave result Wk vanishes at the center of the paramagnetic Bril- 
louin zone. The corresponding Goldstone mode is associated 
with small fluctuations of the ordered moment within the cy- 
cloidal plane. In helimagnets, the spectrum often exhibits a 
second Goldstone mode at the ordering wave vector. This gap- 
less mode is due to fluctuations of the plane of the cycloid. In 
the case at hand, the easy-plane anisotropy generated by the 
DM term forces the cycloidal structure to lie in the be plane 
and creates an excitation gap at the ordering wave vector Q. 

The renormalization of the magnon dispersion within the 
framework of the 1/5 expansion is obtained from the poles 
of the Green's function (Eq. 140^ In Fig.|4]we compare the 
results of the 1/5 expansion with the linear spin wave the- 
ory. It is customary to quantify the effects of the "quantum" 
renormalization of the magnon dispersion by parametrizing 
the latter in terms of "effective" exchange constants J, J' , D 
and comparing them with the "bare" parameters J, J' and D. 

Experimentally, the quantum renormalization is found to 
be rather large, namely j = 1.63(5) and y- = 0.84(9). The 
renormalization of D was not established. The 1 / S expansion 
yields the significantly smaller renormalizations 4 = 1.131, 

= 0.648 and § = 0.72. The difference between the theo- 
retical and experimental values indicates that the leading order 
in a 1/5 expansion underestimates fluctuation effects. On the 
other hand one should note that the 1/5 expansion gives a re- 
sult of 0.031 for the incommensuration, which is very close to 
the experimentally observed value of e = 0.030(2). 

Before turning to a comparison of our results for the dy- 
namical structure factor with the experimental results, we 
briefly review some facts about excitations in helimagnets. 
Generally it is useful to distinguish between three spin wave 
modes. In the case at hand, the "principal" mode = d>k 
is polarized along the a axis and is probed by the 5£ a w com- 
ponent of the dynamical structure factor (Eq. I29> . The two 



"secondary" spin wave modes uj£ = d>k±Q are images of the 
main mode but their momenta are shifted by ±Q. They are 
polarized in the be plane (Eqs. |30] 13 1> . In linear spin wave 
theory, the three spin wave modes give rise to sharp 5 func- 
tions and carry a large part of the spectral weight. 

In addition to the single magnon modes there are multi 
magnon scattering continua. Whenever the magnon disper- 
sion lies within a scattering continuum, the single magnon ex- 
citation gets broadened and acquires a finite line width. On 
the other hand, when the magnon dispersion lies at the thresh- 
old of a scattering continuum, there is no significant decay and 
the single magnon mode remains sharp. 

The unpolarized dynamical structure factor (where the var- 
ious polarizations are added according to Eq.l39t is shown in 
Fig. [3] as a function of energy and momentum for a particu- 
lar "cut" of momentum transfers. The cut along the b direc- 
tion, i.e. from (000) to (010), shows large modulations of 
the dispersion relation due to the strong intra-chain correla- 
tions. Near the ordering wave vector Q the scattering inten- 
sity increases sharply. For momentum transfers perpendicu- 
lar to the chains, (i.e. along the (OI77) direction), the single 
particle modes are seen to be resolution limited. The two 
in-plane modes become degenerate and their dispersions are 
nearly featureless, whereas the out-of-plane fluctuations dip 
to zero energy at (011), in accordance with Goldstone's the- 
orem. Along the (07777) direction the spectrum is symmetric 
across the Brillouin zone boundary. Additional structures due 
to two magnon scattering are clearly visible at higher energies 
along the (O77O) and (O7777) directions. 

In order to illustrate how the spectral weights associated 
with the single-particle excitations are affected by the magnon 
interactions, we have estimated their contributions for each 
polarization to the integrated spectral weights. The total inte- 
grated intensity of each polarization is given by "equal-time" 
correlation functions, 

The one-magnon contribution to the integrated intensity of 
each polarization is then determined by integrating the respec- 
tive correlation function in the vicinity of the single particle 
dispersions. In practice we find that integrating the peaks as- 
suming a Lorentzian form is a poor prescription for strongly 
damped peaks. Instead we numerically integrate the intensity 
over an energy window of three times the width at half maxi- 
mum 

1 /•tnJ+l.SAwJ j^xx 

4 «/u>°-1.5Au,° T 

K£ = t± / duj *±2*L. (42) 

4 A,± ±Q -1.5Au,± ±Q 7T 

The results are shown in Fig. [5] We see that the integrated 
spectral weight is concentrated in the vicinities of the order- 
ing wave vector Q and (0||) and is largely suppressed near 
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Reduced Wavevector If 

FIG. 3: (Color Online) Density plot of the scattering cross-section (Eq. l39> as a function of energy and wave vector. The light (dark) regions 
represent regions of small (large) scattering intensity. The results have been convolved with the experimental energy resolution of the detectors 
(The full width at half maximum is AE — 0.016 meV). The magnetic form factor of copper in Eq. 1391 shows very weak wave vector 
dependence in the regime of interest and therefore was taken to be unity. 32 The filled circles along the (010) direction are the experimental 
position of the most intense peaks in the line shapes taken in the spiral phase (T < 0.1K)i- 



(0r?0) (01)7) (Orw) 





(000) 0.5 (010) 0.5 (011) 0.5 (000) 

Reduced Wavevector r\ 

FIG. 4: (Color Online) The renormalization of the spin wave spec- 
trum. The solid and dashed lines are results of the 1/S expansion iDk 
(Eq, l40> and the linear spin wave dispersion lo^ (Eq. ll6t . The cut in 
the paramagnetic Brillouin zone runs from (000) to (010) to (011) 
and back to the center of the Brillouin zone (000). 



FIG. 5: (Color Online) Upper panel (A): The integrated spectral 
weights for the three polarizations (Eq. 14 It as functions of momen- 
tum transfer. The principal lu° and secondary uj + , uj~ spin wave 
modes are defined in the text. Lower panel (B): The ratios of the 
spectral weights of the single-particle excitations of each polariza- 
tion to their respective integrated intensities (Eq, l42t . 



the r point. The weights associated with single magnon exci- 
tations are strongly suppressed for the secondary modes. This 
is a consequence of the non-collinearity of the magnetic or- 
der. The in-plane modes are significantly damped as a result 
of the coupling between longitudinal and transverse fluctu- 
ations. Such a coupling is not present for the out-of-plane 
mode and therefore the principal mode carries generally more 
spectral weight. Nevertheless the fraction of spectral weight 
associated with single-particle excitations decreases signif- 
icantly whenever the renormalized spin wave dispersion is 



pushed upwards in energy for a given momentum. For in- 
stance, near wave vector (0, 0.8, 0) the principal spin wave 
mode lies within the two magnon continuum and as a result 
less than 50 % of the spectral weight is attributed to the one- 
magnon excitation. 

The scattering intensity can also be studied by performing 
a wave vector average, 

k £i 
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By the frequency sum rule, the scattering intensity (Eq.I43> in- 
tegrated over all energies (including the elastic Bragg peaks at 
u> = 0) has to equal 5(5+1). However, this sum rule does not 
hold "order by order" in perturbation theory. For instance, the 
total intensity within linear spin wave theory exceeds the sum 
rule by A5(l + 2A5). Bearing this caveat in mind, the sum 
rule is a useful tool for comparing the one and two magnon 
contributions as well as analyzing the shift in spectral weight. 
In Fig.|V] we plot the scattering intensities as functions of en- 
ergy within linear spin wave theory and the 1 / S expansion. In 
linear spin wave theory the integrated intensity exhibits cusps, 
which are associated with van-Hove singularities in the sin- 
gle particle density of states. In the 1/5 expansion such sharp 
features are absent. Above approximately 0.5 meV the one 
magnon contribution vanishes and the scattering intensity is 
entirely due to multi magnon states. 

To quantify the shift of the spectral weight we calculate the 
first moment of the normalized scattering intensity (ui). We 
find that the linear spin wave theory value (ui) = 0.35 meV 
is renormalized upwards to (ui) = 0.40 meV in the 1/5 ex- 
pansion. This observation is in line with the expectation that 
the higher orders of the 1/5 expansion induce a transfer of 
spectral weight to higher energies via multi magnon scatter- 
ing processes. In fact, as shown in Fig. [V] the two magnon 
contribution to the overall intensity is 29% in linear spin wave 
theory but 46% in the 1 / S expansion. 
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FIG. 6: (Color Online) The scattering intensity as a function of en- 
ergy (Eq. l43t for LSW theory (dashed line) and LSW+l/S expansion 
(solid line). The contributions of the scattering continua is shown 
using thin lines. 



A. Excitation line shapes 

In order to exhibit the properties of the dynamical structure 
factor in greater detail we have generated a series of scans in 
k uj space. The inelastic neutron scattering measurements on 
CS2CUCI4 were not performed at constant momentum transfer 
but followed various trajectories in energy-wave vector space. 
We have generated our theoretical scans using the known pa- 
rameterizations of the scans A to J of Ref. 4 in k — u> space, 
which we summarize in Table. H] We refer the reader to Ref.|4j 



FIG. 7: (Color Online) Density plot of the scattering cross-section as 
a function of energy and wave vector along the cut (Ofc§). The light 
(dark) regions represent regions of small (large) scattering intensity. 



for further details. The various scans are shown in Fig. [8] 
Also shown are the regions in which significant magnetic scat- 
tering is observed experimentally and the location of the main 
peaks. For comparison we plot the principal and secondary 
spin wave dispersions obtained from the 1/5 expansion. As 
we have already emphasized, the 1/5 expansion underesti- 
mates the quantum renormalization of the exchange constants 
and as a result the agreement of the calculated spin wave dis- 
persions with the main peaks observed experimentally is poor. 

The experimental energy and momentum resolutions have 
been accounted for to make contact with experiment. We 
find that the effects of the finite energy resolution of AE — 
0.016 meV are generally outweighed by the effects of the fi- 
nite momentum resolution. This is a consequence of the large 
modulation of the spin wave dispersion along the chain di- 
rection, i.e (OfcO), (whose slopes can reach ~ 1.6 me V), 
which causes an amplification of the effects of the momentum 
resolution. Given that the spin waves are nearly dispersionless 
along the (001) direction, we have only taken into account the 
spatial resolution along the chain direction. 

To illustrate this point, let us consider the results for scans 
B, E, G and H shown in Fig. [9] The insets of panel (4) show 
the results of both linear spin wave theory and the 1/5 expan- 
sion for a hypothetical energy resolution of AE — 0.002 me V 
which has been introduced to make the various delta func- 
tion peaks visible (the momentum resolution is set to zero 



# fc(rlu) Krlu) ka_ 

A —0.389 + 0.189-ZJ — 0.016-E 2 0.05 

B -0.30 + 0.189S - 0.015S 2 0.02 

C 0.21 + 0.297S - 0.026S 2 0.25 

D 2.11 + 0.29E - 0.025S 2 0.95 

E -0.33 + 0.19S-0.015£ 2 0.78 + 0.37E - 0.03E 2 1 
F -0.39 + 0.19S - 0.02S 2 1.66 + 0.37E - 0.035S 2 1 
G 0.5 1.53 - 0.32E - 0.1E 2 1 

H 0.28 + 0.29E - 0.025E 2 1.205 1 

J 0.47 1.0 - 0.45£ 1 

K 0.29 + 0.29E - 0.03£ 2 0.77 - O.UE + 0.013S 2 1 



TABLE I: Parameterization of energy-momentum scans performed 
in Ref.0: the momentum transfers k = (h,k,l) are parameterized 
in terms of the energy transfer E (in meV). k a is a measure of the 
polarization factor. Given that the weak interlayer coupling is ne- 
glected, h is not needed for the purpose of our calculation. 



Afc = 0). First we consider the results for scan H (Panel 4 
of Fig. [9}. Linear spin wave theory predicts peaks at approxi- 
mately 0.27 meV and 0.37 meV corresponding to the degen- 
erate spin wave modes and to u~ respectively. The 
1/5 correction yields a slight upward shift in the energy ol 
these peaks. In both linear spin-wave and 1/5 calculations, 
the two magnon scattering continuum is found to carry nearly 
a quarter of the integrated spectral weight. Taking into ac- 
count the finite momentum resolution (the width at half max- 
imum is Afc = 0.057) we find that the sharp peaks get broad- 
ened very significantly as is shown in panel (4). The dynam- 
ical structure factor now exhibits an extended continuum in 
which the single-particle excitation can no longer be resolved 
and merges smoothly with the two magnon continuum. This 
result is qualitatively similar to the experimental observations 
shown for comparison in panel (2) of Fig. [TO] 

Next we turn to scan G (Panel 3 of Fig. |9}, which probes 
the vicinity of the wave vector (0,0.5,1.5). Experimen- 
tally a resolution-limited peak is observed at an energy of 
0.107(10) meV in this region of intense scattering, see panel 
( 1 ) of Fig El| However, about two thirds of the spectral weight 
is associated with a scattering continuum at higher energies. 
Both linear spin wave theory and the 1/5* expansion predict 
sharp peaks in this region of the Brillouin zone. The 1/5 ex- 
pansion gives a spin wave peak at u° = 0.18 meV carrying 
nearly half of the spectral weight and two further peaks at en- 
ergies around 0.25 meV corresponding to the two secondary 
spin wave modes. The two magnon scattering continuum ex- 
tends up to 0.9 meV and carries nearly a quarter of the spectral 
weight. We emphasize that, in contrast to u> , the principal 
mode cu° is close to a saddle point and therefore is nearly dis- 
persionless. In Panel 3 the finite energy and momentum reso- 
lutions are taken into account. We see that the almost disper- 
sionless main mode remains sharp but the secondary modes 
can no longer be resolved and are found to merge with the 
two magnon continuum. Irrespective of the discrepancies be- 
tween the results of the 1/5 expansions and the experimental 
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FIG. 8: (Color Online) The dispersion relation of magnetic excita- 
tions. The shaded regions labeled with capital letters A through K 
indicate scan directions (the line thickness indicates the wave vector 
averaging). The filled symbols are the main peaks in the line shape as 
determined experimentally in the ordered phase (from Ref.0). The 
dotted area indicates the extent of the scattering continuum. The 
open circles and squares are respectively the upper and lower bound- 
ary of the scattering continuum as determined experimentally. The 
upper thick dashed line is a guide to the eye. The thin solid line is 
the experimental fit to the principal mode using effective parameters 
(J = 0.61 meV, J' = 0.107 meV). The thick solid, dashed and dash 
dotted lines are respectively the 1/S results for the principal (w£) 
and secondary uj^) modes determined from Eg. 1401 



data, our calculation suggests that the lower boundary of of 
the measured scattering continuum in scan G could be due 
to unresolved transverse magnons. Such a scenario had been 
previously considered and ruled out on the basis of the small- 
ness of the ratio I se c/Ipri of spectral weights of the secondary 
modes to the principal mode predicted by linear spin wave 
theory. 4 However, the results of the 1/5 expansion show that 
spin wave interactions lead to an enhancement of this ratio for 
the G scan. 

Next, we examine scan E (Panel 2), which probes wave vec- 
tors near k = (0, —0.25, 1). Linear spin wave theory predicts 
coherent peaks at w° = 0.35 meV for the principal mode and 
at u~ = 0.44 meV and oj + = 0.33 meV for the secondary 
modes (see the inset in Panel 2). The two magnon scatter- 
ing continuum is relatively weak and carries only about 23 
% of the total spectral weight. In the framework of the 1/5 
expansion the principal mode is pushed upwards in energy to 
uj° = 0.42 meV and occurs very close to the secondary mode 
oj~ = 0.45 meV. The other secondary mode uj + is shifted 
very significantly to 0.39 meV, but carries only a minute frac- 
tion of the spectral weight. The two magnon continuum is also 
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FIG. 9: (Color Online) Scattering cross section. The numbered panels (1-4) correspond to energy scans B, E, G and H respectively. The data 
has been convolved with the energy and spatial resolution. AE = 0.016 meV for all plots and Afc = {0.035, 0.039, 0.085, 0.056} for plots 
(1-4). The insets show the results of linear spin wave (LSW) theory and the 1/5 expansion (LSW+l/S) for Afc = 0, AE = 0.002 meV. 
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FIG. 10: Observed neutron scattering lineshape in scan G (1) and 
H (2) (data from Ref.0) in the ordered phase (T < 0.1K). In scan 
G, the shaded area represents the the 1/S calculation (Eq. (39)) con- 
volved with the experimental resolution. 



shifted upwards in energy and carries approximately a quarter 
of the total spectral weight. Once again the spin wave disper- 
sion is close to a saddle point and as a result the effects of the 
finite momentum resolution are small. The main feature in the 
structure factor is a broad peak formed by the two unresolved 
uj~ and ui° modes. This is quite similar to what is observed 
experimentally (Fig. 5(E) of Ref. |4J). It is then tempting to 



speculate that the experimentally observed single peak is a re- 
sult of the accidental near degeneracy of the lo~ and w° modes 
in the vicinity of k = (0, —0.25, 1). This would explain both 
the absence of the cu~ peak in the experimental data and the 
anomalously large intensity of the observed peak. 

In Panel 1 of Fig. [5] we plot the dynamical structure factor 
for scan B near the wave vector (2, —0.25, 0). Here the po- 
larization factor (k a ) in fl39i leads to a strong suppression of 
the out-of -plane fluctuations and the scattering is almost en- 
tirely due to the in-plane lo spin wave modes. The magnon 
interactions renormalize lo + upwards in energy to approxi- 
mately 0.42 meV, whereas the uj~ mode disappears in the two 
magnon scattering continuum. A careful analysis shows that 
the narrow peak at 0.55 meV is not due to a single-particle 
excitation but is a feature in the two magnon scattering con- 
tinuum. 

The dominant contribution to the dynamical structure fac- 
tor in scan A in the vicinity of the wave vector (1.5, —0.3, 0) 
comes from in-plane fluctuations because the polarization fac- 
tor k a suppresses out-of plane fluctuations. As can be seen 
in Fig. the magnon interactions lead to a spectral weight 
transfer to higher energies. The peaks near 0.8 meV and 
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FIG. 1 1 : (Color Online) Calculated scattering cross sections in linear 
spin wave theory (LSW) and 1/S expansion (LSW+l/S) for scans 
A and C. Panels (2) and (4) show the LSW+l/S results with instru- 
mental resolutions of AE = 0.016 meV and Afc = 0.02 (for scan 
A) and Afc = 0.04 (for scan C) taken into account. See Fig. 5 from 
Ref.0to compare to experimental data. 

0.85 meV can be traced back to single-particle poles in the 
Green's function. These poles are unphysical and are a result 
of the uncontrolled nature of the 1 / S expansion for small val- 
ues of 5. It is easily seen from the Dyson equation J2 It that 
a large self-energy at a given wave vector can lead to "extra" 
poles in the Green's function at high energies above the two 
magnon continuum. The inclusion of higher order terms in the 
1/5 expansion would provide decay mechanisms at all ener- 
gies and lead to a broadening of these high-energy peaks in 
the dynamical structure factor. 

Last but not least let us consider the vicinity of (0.8, 0.4, 0) 
(scan C). As is shown in Fig. the principal spin wave 
mode lu° is renormalized down to a slightly lower energy 
of approximately 0.42 meV. The lu + mode, which occurs at 
0.35 meV in linear spin wave theory, disappears entirely in 
the two magnon-continuum. The feature near 0.60 meV can 
again be understood in terms of an enhancement of the two 
magnon density of states. Comparing with the neutron scat- 
tering data (Fig. 5(C) of Ref. 4), the structure factor shows 
features quite similar to the experimentally observed contin- 
uum. However, the scattering continuum occurs at energies 
nearly 0.10 meV lower than what is observed experimentally. 

VI. CONCLUSIONS 



value and geometrical frustration. 

We found that the results of nonlinear spin wave theory ex- 
plain on a qualitative level many of the features observed in 
neutron scattering experiments. We find a strong scattering 
continuum in the dynamical structure factor similar to the ex- 
perimental observations. Our calculations suggest the pos- 
sibility that some of the spectral weight at the low-energy 
boundary of the experimentally observed scattering contin- 
uum in scan G could be due to single particle excitations that 
are unresolved. 

In the vicinity of saddle points of the spin wave disper- 
sion relation the single-particle excitations are only weakly 
affected by the instrumental resolution and hence exhibit 
sharper peaks in the dynamical structure factor. 

In spite of the qualitative agreement of the theory with 
experiments, crucial discrepancies remain. First and fore- 
most, nonlinear spin wave theory fails to account for the large 
"quantum renormalization" of the main exchange parameter. 
This indicates that (to order O(S )) the 1/S* expansion still 
underestimates the effects of quantum fluctuations. Further- 
more, there are significant quantitative differences between 
our calculations and the experimentally observed structure 
factor. One may speculate that a better agreement with ex- 
periment could be achieved by taking higher-order terms in 
the 1/5 expansion into account. 

The main lesson to be learned from our calculations is that 
CS2CUCI4 falls somewhere in between the two theoretical 
scenarios that have been proposed previously. Our analysis 
shows that the physics of order plays an essential part in un- 
derstanding the dynamic response CS2CUCI4 at low tempera- 
tures: a large fraction of the spectral weight is carried by spin 
wave modes, which occur over a large range of frequencies. 
This is a strong indication that a putative spin-liquid ground 
state is plainly not a good starting point for the description of 
the ordered phase of CS2CUCI4. On the other hand we have 
seen that (in low orders in 1/5) nonlinear spin wave theory 
significantly underestimates the effects of quantum fluctua- 
tions and hence expansions around the ordered state also fail 
to account for the experimental observations. 

Nonlinear spin wave theory can also be applied to inves- 
tigate the effects of magnetic fields. It is known that in the 
presence of a field linear spin wave theory is generally a very 
poor approximation as it excludes all-important magnon de- 
cay processes^ A self-consistent study of magnetic field ef- 
fects in CS2CUCI4 is currently under way. 34 During comple- 
tion of this work, we became aware of a parallel effort which 
reaches similar conclusions^ 



In this work we have used nonlinear spin wave theory to de- 
termine the dynamical structure factor in the ordered phase of 
the spin-1/2 helimagnet CS2G1CI4. We have taken into ac- 
count the first subleading contribution in a 1/5 expansion, 
which incorporates interactions between magnons and gen- 
erates magnon decay processes as well as multi magnon scat- 
tering continua. Both effects are particularly pronounced in 
CS2CUCI4 due to the non-collinear spin ordering, the low spin 
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